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We study non-perturbatively, via the Schwinger-Dyson equations, the leading 
infrared behavior of the pressure in the ladder approximation. This problem is dis- 
cussed firstly in the context of a thermal scalar field theory, and the analysis is 
then extended to the Yang-Mills theory at high temperatures. Using the Feynman 
gauge, we find a system of two coupled integral equations for the gluon and ghost 
self-energies, which is solved analytically. The solutions of these equations show that 
the contributions to the pressure, when calculated in the ladder approximation, are 
finite in the infrared domain. 
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I. INTRODUCTION 



Relativistic field theories at finite temperature have been actively studied in the past 
years, because of their relevance to the theory of the early universe and to the quark-gluon 
plasma which may be created in heavy ion collisions [[[]]. The study of such physical problems 
requires the computation of the thermodynamic potential, from which other thermodynamic 
properties, like the pressure, may be determined. It was pointed out by Linde 0] that the 
thermodynamic potential of the Yang-Mills theory cannot be calculated beyond the fifth 
order of the coupling constant. This is related with the infrared singularities of the non- 
abelian gauge theories at finite temperatures, which arise from the fact that the magnetic 
mass vanishes at least up to second order in perturbation theory. The infrared problem in 
these theories at finite temperature is qualitatively different from that at T = 0, manifesting 
itself by the presence of infrared power divergences in higher orders of the perturbative 
expansion. 

One may hope that the use of non-perturbative methods can throw some light on this 
problem. Such methods generally start with a discussion of the relevant Schwinger- Dyson 
equations ||. It is well known that this set of integral equations is not closed since an n- 
point function is generally related to the n + 1 and n + 2 functions. This is what makes the 
solutions of such problems exceedingly difficult. For this reason one must resort in practice 
to some approximations, which result from truncating in some way the set of the Schwinger- 
Dyson equations. This non-perturbative approach has been used by Jackiw and Templeton 

to show how super- renormalizable interactions might cure their infrared divergences. It 
has also been employed by Mandelstam || to analyze the infrared behavior of the gluon 
propagator, using the ladder approximation at T — 0. He showed that the ensuing set of the 
Schwinger-Dyson equations does provide confinement in the strong coupling regime. 

The purpose of this work is to study, in the ladder approximation, the leading infrared 
behavior of the pressure in the Yang-Mills theory at high temperatures. Our approach, via 
the Schwinger-Dyson equations, is similar in spirit with that advocated by Kajantie and 
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Kapusta 0. Namely, we regard the simplified set of these equations just as a convenient 
procedure for summing particular classes of an infinite number of diagrams in the weak 
coupling regime. In order to regularize the theory in the event that no magnetic masses 
are generated, we shall put an infrared cut-off A on the momenta integrations which arise 
when calculating the contributions to the pressure. Our aim is to study the possibility of 
cancellation of the infrared singularities in the limit A — > 0, when summing the whole class 
of diagrams which contribute in the ladder approximation. 

As we shall see, even this more modest task is non-trivial and we begin by considering 
in Section |IJ the scalar g<p 3 theory in six dimensions. This model has some similarities with 
the Yang-Mills theory, such as a dimensionless coupling constant and asymptotic freedom. 
Strictly speaking, in this case one could employ the re-summation method developed by 
Braaten and Pisarski [[/J , because a scalar thermal mass is generated in lowest order. Indeed, 
such a procedure leads to an infrared finite expression for the pressure || . However, we will 
use here instead the method described above, in order to illustrate in a simple way some 
relevant features which will appear in the Yang-Mills theory. We consider the truncated set 
of Schwinger-Dyson equations which yield consistently the scalar self-energy in the ladder 
approximation. The ensuing integral equation is equivalent, under certain conditions to be 
discussed later, to a second order differential equation which can be solved analytically in 
terms of modified Bessel functions. Using this solution, we show that the corresponding 
contribution to the pressure obtained in this approximation, is finite in the limit A — > 0. 



In Section |TTTj we consider the leading infrared behavior of the SU (N) Yang-Mills the- 
ory at high temperatures , in the ladder approximation. We work in the Feynman gauge 
and obtain a system of two coupled Schwinger-Dyson integral equations, linking the gluon 
polarization tensor and the ghost self-energy function. Using the same reasoning as in the 
previous section, we show that this system is equivalent to a fourth order differential equa- 
tion which describes in this approximation the leading infrared behavior. The solution of 
this equation is expressed in terms of modified Bessel functions with complex argument and 
the corresponding complex conjugate functions. We find that the relevant contributions to 



the pressure, as calculated in the ladder approximation, are finite in the infrared limit A — > 0. 
Some of the mathematical details which arise during these calculations are given in the Ap- 
pendices A and B. Finally, in Appendix C we consider the possibility of the cancellation of 
the infrared divergences in the pressure, beyond the ladder approximation. 



II. THE SIX-DIMENSIONAL SCALAR THEORY 

We consider here the scalar theory with a g<p 3 interaction at the temperature T. In 
this theory the formula giving the pressure in the ladder approximation is rather simple 
and can be expressed directly in terms of the self-energy function EL The Schwinger-Dyson 
equation for the scalar self-energy in the ladder approximation is shown in Fig. (1). Note 
that the diagrams contributing to II have the same combinatorial factor ~ to all orders in 
perturbation theory. The corresponding contributions to the pressure are represented in Fig. 
(2.a). 

In order to derive the expression giving the pressure, we use the relation between its 
functional derivative and the self-energy function [|I|]: 

where T>q is the free particle propagator and 1PI indicates that only one particles irreducible 
diagrams contribute to EL With the help of the graphical representation given in Fig. (2. a), 
it is easy to see that (|l|) implies the relation: 

1 ~ 



(2) 



The factor | arises because in general there are only four ways we can cut the diagram 
(2. a) in order to get IT in the ladder approximation. The only exception occurs in the lowest 
order, where the factor is | instead. The graphical representation of this formula is shown 
in Fig. (2.b). 

The sum over po should be taken over even frequencies po = 2ixnT . However, the 
dominant infrared contributions arise only from the terms with zero frequency (n = 0). 



When calculating the contributions from eq. (0) corresponding to this mode, we can put 
an ultraviolet cut-off of order T on the momentum integration. This cut-off arise naturally 
when summing over all modes n. As mentioned in the Introduction, we shall also put a 
cut-off A in order to regularize the perturbative infrared divergences. In this way we find 
that the leading infrared contributions to the pressure are given by: 

T rT d 5 p 1 ~ 



P = Pn 



r n(p,p = o) + 



(3) 



/a (2tt)V 

where dots denote additional subleading and infrared convergent contributions. 

Then, using the Schwinger-Dyson equation for the scalar self-energy function we arrive 
at the following integral equation for H(p,p = 0): 



n(p) = u(p)+ P 2 



9 2 T 



d 5 k 



i n(fc) 

2k 2 k A 



(4) 



(2tt) 5 A (p + k) 2 

In order to solve this integral equation analytically, we will try to convert it into a 
Volterra type. This is possible provided we adopt the following maximization procedure. 
Let us define: 



(P + k) 



2 

max 



p 2 for p > k 
k 2 for k > p 

Using this approach, we then obtain from eq. (01): 



(5) 



n ma x(p) = -p 2 + 



12tt 3 



p T 1A 3 1 fP , 1TT .,. f 1 dk^ .,. 
3 2 6 p 2 ja J P k 



T dk. 



Comparing this with the original equation (|J), we note that the first iteration of 
A — » yields: 



2 ^ ff^ 

^ +24^ 



2,2 | lii 



3T 



(6) 
with 

(7) 



a result which completely agrees in the high temperature limit with the leading contribution 
which follows from (0J). To next order, our procedure yields in this regime results which 
are in satisfactory agreement with the ones obtained from equation (HJ), which becomes 
increasingly cumbersome to handle . Since we wish to study only the main features of the 



leading infrared behavior, it will be sufficient for our purpose to restrict our attention to the 
Volterra equation (|6|). For simplicity of notation we shall drop in what follows the suffix max 
appearing in (|6]). After successive iterations, we encounter the presence of infrared divergent 
terms in the perturbative series of II: 

/ r f~* \ / T \ ^ / \ T / T \ ^ 

„2 i a rp2 „2rp2 / 1 \ i „3rp2 1 \ i„(P\ ' „ 4^2 I 1 \ P 



m = V + _r - *T> - + oft- - In ^ - -a>T< - f + 



2 

where a = t?t 



12tt 3 • 

With the help of this result and using equation (|3]), we find that the perturbative expan 
sion of the pressure exhibits power infrared divergences given by: 



div 



A 6 \XJ 36 \\J 540 \XJ 



(9) 



967T 3 

Despite the fact that equation (^|) cannot be solved perturbatively because the iterations 
yield infrared divergent terms, we will show that nevertheless, a non-perturbative solution 
does exist. To this end, it is convenient to define: 

olT „. , n . 

x = — ■ f(x) = - 10 
p jr 

Then (||) reduces, dropping the inhomogeneous A 3 term, to the following integral equa- 
tion: 



f(x) = -1 + X - - f - x 2 f %f{y) - * 4 *Lf(y) (11) 
3 2a Ja y z Jx 



Perturbation theory now corresponds to solving (|TT1) by a power series in x, a procedure 
which, as we have seen, yields power infrared divergences. 
But the differential equation which follows from (|TT|): 



x 2 f" - bxf + (8 - 2x)f = x - 8 (12) 

has a well behaved solution. The complementary one involves modified Bessel functions || 
and two constants: 
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f c (x) = Ax 3 I 2 {VSx) + Bx 3 K 2 (V8x) 



(13) 



Using standard methods we can then find a particular solution in terms of these modified 
Bessel functions. The general solution, consisting of the sum of f c (x) and the particular one 
has the form: 



/(*) 



A + 2 



+ 



B 













1 







X 3 Io 



x 3 K, 



lx) + 



ix) 



(14) 



Substituting this into the integral equation (|TTD , fixes the constants A and B in terms 
of the parameters ^ and a. Since I 2 grows exponentially at large x, it would produce 
a divergence in the integral equation as A —>■ 0, unless A = in this limit. Indeed, the 
consistency conditions on f(x) at 1 = y and x = a, demand A to vanish when A — > 0, 
while B becomes: 

1 



B 



2 i 3 (^)/;^g-i)^)-jA 



(15) 



Then, since K2 decreases exponentially at large x, it is not difficult to show that f(x) 
is a decreasing function of x in this domain. Because large values of x correspond to small 
values of the momenta [see eq.([T0|)1, this behavior leads to a convergent integral in ([|) as 
A — ► 0. Thus, the pressure calculated in the ladder approximation remains finite in the 
infrared domain. This also happens in the case of the Yang-Mills theory to which we now 
turn. 



III. THE YANG-MILLS THEORY 

We now consider the leading infrared behavior of the pressure in thermal Yang-Mills 
theory. Since the IIoo(p ~^ 0,po = 0) component of the polarization tensor is non-zero, 
the longitudinal gluons are screened at large distances. As a result, the infrared problem is 
related in this case only to the behavior of the transverse part of the gluon propagator. Using 
the fact that at finite temperatures, the leading infrared contributions arise from terms with 
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zero frequencies, we can reduce our problem to a study in 3-dimensional Euclidian Yang-Mills 
theory. 

Although the analysis is now more complicated due to the presence of the 4-gluon cou- 
plings and ghosts, it proceeds in parallel with that described previously in the Section [IT]. In 
particular, in order to derive the expression for the pressure in the ladder approximation, we 
start from a basic relation which is the analogue of eq. (p]). By a reasoning similar to that 
used in deriving eq. (ff) , we find that the expression giving the leading infrared contributions 
to the pressure has the form: 

p = p » " T [ (84 \^ tip) ~ ^ {p) + 5iJ«* + ' ' ] (16) 

This relation is represented graphically in Fig. (3). Here IT and £ denote respectively the 
gluon and ghost self-energy functions, which must be calculated consistently in the ladder 
approximation. The vertex V stands for the bare gluon four-point function and dots indicate 
infrared convergent contributions associated with two-loop diagrams. Note that expression 
(PS|) is much simpler than the exact formula giving the pressure in the Yang-Mills theory 

We now proceed to investigate, in the Feynman gauge, the infrared behavior of the gluon 
polarization tensor and the ghost self-energy function. The corresponding Schwinger- Dyson 
(S-D) equations in the ladder approximation are shown in Fig. (4). 

It is well known that any approximation of the S-D equations may impose severe con- 
straints in the case of a gauge theory. The reason is that the Ward identities, which reflect 
its underlying gauge invariance, are satisfied only when we take into account all the rele- 
vant contributions, order by order in perturbation theory. In particular, the transversality 
property of the polarization tensor is guaranteed only by the full set of S-D equations. Since 
it is impossible to implement this program in practice, we follow the procedure adopted by 
Mandelstam ||, neglecting the longitudinal terms which might arise in connection with the 
approximate set of S-D equations. This is obviously correct to lowest order in perturbation 
theory, where our polarization tensor is manifestly transverse. In higher orders, some of 
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the contributions associated with the exact 3-point and 4-point gluon vertices will cancel 
these longitudinal terms. Since in the ladder approximation the vertex corrections are also 
neglected, the above procedure is justified and, as we shall see, a consistent solution of the 
approximate set of S-D equations does exist. 

Thus, using the transversality property of the gluon polarization tensor 11°* = 5 fly, we 
can write it in the form: 

%= (^-^r)(nT + P 2 ) (17) 

In this way, we arrive at the following set of coupled integral equations which relate the 
gluon polarization tensor to the ghost self-energy S afe = <5 afc (£ + p 2 ): 

_ 2 _llg 2 NT rd 3 k g 2 NT r d?k 2(2k 2 + 3p 2 )(p • k) 2 - p 2 k 2 {6p 2 + 7k 2 ) 

T ~ ~ P 48tt 3 J ~W ~ 16tt 3 J ~¥~ p 2 (p + k) 2 

g 2 NT fd 3 k TT n , g 2 NT f d 3 k 5p 2 k 2 - 3(p • k) 2 , 2 )2ntt /7 . 

I -—Tl T (k) + g— - — F v f ' (p 2 + k 2 )U T (k) + 



24tt 3 J k A iy ' 8tt 3 J k 6 p 2 (p + k) 2 
g 2 NT r d 3 k p 2 k 2 - (p ■ k) 



Stt 3 



k A p 2 (p + k) : 



-E(fc) (18) 



2 ^ 2 iVT /• rf 3 fc (k-p) 2 /iVT r rf 3 fc (p-k) 
S (p) = -p - — — / — — + — — / — S(fc) - 



8tt 3 J k A (p + k) 2 8tt 3 J k 4 (p + k) 2 
^ 2 A^T /■ d 3 k p 2 k 2 - (p ■ k) 2 



n r (fc) (19) 



8tt 3 7 fc 6 (p + k) 2 
The first iteration of this system yields, as required, the lowest order perturbative con- 
tributions to the gluon and ghost self-energy functions. In higher orders the perturbative 



series of these self-energy functions lead, via equation (|1^), to the presence of power infrared 
divergences in the perturbative expansion of the pressure. These features are rather similar 
to the ones exhibited in Section || by the thermal scalar field [see eq. @ and (0)]. 

Analogously to the scalar case, we will now show that a well behaved non-perturbative 
solution for the above set of integral equations does exist. In order to be able to solve it 
analytically, we will use consistently in the numerators and denominators of equations (|IBD 
and (0), the maximization procedure described in the previous case. As we have seen, in 
the high temperature domain this procedure simplifies the equations, without modifying the 
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qualitative features of their solutions. Then, it is convenient to define the dimensionless 
quantities: 



- 9 2N - aT fort 

a = —— ■ x = — 20 

67r z p 



F(x) = ^- ■ G(x) = - 2 (21) 
p z p z 

In this way, after performing the angular integrations, we find the following set of Volterra 

integral equations: 



F(x) = -1 - ^x - Ex 2 + x 2 f %G{y) + x 4 %G{y) (22) 
i/ Ja y Jx y 











+ x- f 




'2) 


J a 



aT t 

— ay 



My) (23) 

x y 



where 



1 11 r^r dy 



H = — + — -^F(y) (24) 
2a 2 Ja y 2 KyJ V ; 

Perturbation theory corresponds to solving ([22|) and (|23"D by a power series in x. On 



the other hand, the relevant momenta in the infrared domain are such that p < aT, which 
correspond to large values of x. 



We now consider the differential equations which follow from (|22|) and (E3 



39 

x 2 F"(x) - hxF\x) + 8F(x) + 2xG{x) = -8 - — x (25) 
xG"{x) - G'{x) - 2F{x) = -1 (26) 

To obtain the solution of this set of coupled equations, we use (|25| ) to express G(x) in 

function of F(x) and its derivatives. Substituting this into fl26|), we get the following fourth 

order differential equation: 

4 12 24 /24 4 \ 2 24 

F""(x) - *F»'(x) + - 2 F"{x) - ^F'(x) + ( q + A ) F (x) = ~ 2 - ~ 4 (27) 

which describes the infrared behavior of our system. 

The complementary solution of the homogeneous part of equation fl2~T|) can be expressed 
in terms of the modified Bessel functions with complex argument. Using the functional 
relations satisfied by these functions §J it is straightforward to show that: 
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F c (x) = ix 2 [C\I 2 (V8ix) + C 2 K 2 (V8ix)} + complex conjugate 



(28) 



where C\ and C 2 are complex constants to be determined. 

In order to obtain a particular solution of (|27|), we employ the method of variation of 
the parameters, which is described in Appendix A. The general solution, consisting of the 
sum of F c (x) and F(x) [equation (|A7|)1 is given by: 



F(x) 



ix 



+ 



Ci + 



lit) 



x / 1 1 2 N 
C 2 + / dt[---)I 2[ 

a 



lit) 



K 2 (V8ix) \ + c.c 



(29) 



J 2 t 4 . 

Since I 2 grows exponentially for large values of its argument, this result would produce 
a divergence as A — > in the integral equations ( P^| ) and (^), unless C\ vanishes. In 
Appendix B, we discuss the consistency conditions which follow from these equations in the 
limit A — > 0. We find that, indeed, these conditions require the vanishing of the constant 
C\. Futhermore, they fix the constant C 2 in terms of a, but the explicit expression has a 
rather complicated form involving the modified Bessel functions [see eq. (|B2|)]. 

From the behavior of F(x) and G(x) for large values of x, given by (|B5| ) and (|B6|) , and 
using eqs. (|20|j2l|) , we see that Ht(p) and S(p) will be proportional to p 2 for small values of 
the momenta: 



n r (p) 



E(p) 



39 



TP 



[p« aT] 



(30) 



2 1 7 yi ' 2A 1 

This behavior will ensure the convergence of the integral in eq. (|16D as A — > 0. Hence we 
find that the leading contributions to the pressure, when calculated in the ladder approxi- 
mation of the SU(N) Yang-Mills theory, are finite in the infrared domain. 

In conclusion, we stress that our analysis was restricted to the ladder approximation, 
in which case the set of Schwinger-Dyson equations remains linear, allowing for an analytic 
solution. However, there are many other important contributions which must be taken 
into account when studying the infrared behavior of the pressure. These include the set 
of diagrams with crossed ghost and gluon lines, which can be reduced in the large A^- 
limit to the class of planar diagrams . Furthermore, one must also consider additional 
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self-energy and vertex corrections, which are relevant for the appearance of the effective 
coupling constant g(T). It is conceivable that the inclusion of these contributions, may 
indicate the presence of new collective phenomena associated with the infrared behavior of 
the pressure. Nevertheless, in view of the above results, we believe that the cancellation of 
the infrared divergences is a possibility which deserves further investigation. 
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APPENDIX A 



Here we derive a particular solution of eq. fl27|), using the method of variation of the 
parameters. The complementary solution was given in eq. (^3), so we assume a particular 
one of the form: 



F(x) = ZAx)x 2 l 2 { 



nx) 



Z 2 (x)x 2 K 2 ( v8ix) + c.c. 



(Al) 



where Z\{x) and Z 2 {x) are two complex functions to be determined consistently. We have 
to our disposal four real functions, but only one condition, so that we are free to impose 
three more conditions. We choose these additional conditions as follows: 



Z[(x)x 2 I 2 (V8ix) + Z' 2 (x)x 2 K 2 (VSix) + c.c. 



Z[(x) 



x 2 h( 



nx 



+ Z'(x) 



x 2 K 2 (V8ix) 



Z[(x) x 2 I 2 (V8ix) + Z' 2 (x) x 2 K 2 (V8ix) 



c.c. 



+ c.c. 









(A2) 
(A3) 
(A4) 



Then, the expressions giving F'(x), F"(x), F"'(x) and F""(x) simplify considerably. 
We now impose the basic condition that ( |A1| ) be a solution of ( p7|) . Thus we substitute 
these expressions into ([27]) and obtain an identity. Since x 2 I 2 (\/8ix), x 2 K 2 (^/8ix) and their 
complex conjugates are solutions of the homogeneous equation, this identity reduces to: 

2x 2 - 24 



Z[(x) 



x 2 h 



nx 



+ Z'(x) 



x 2 K 2 ( 



■nx 



+ c.c. 



(A5) 



Together with the auxiliary conditions ( |A2|) , ( |A3| ) and (|A4f ), we have a system of equa- 
tions which determinate the complex functions Z[(x) and Z' 2 (x). To this end we use the 
functional relations satisfied by K 2 (z) and I 2 (z), together with the fact that their Wronskian 
is -. Then, after a straightforward calculation we obtain: 



Z[(x) = -i[- 2 - l ^)K 2 {^%ix) 



Z' 2 {x) 



^-•2 J 



or 1 / 



nx) 



(A6) 



Integrating these relations and substituting the result into (|AT|), we find that the partic- 
ular solution F(x) can be expressed in the form: 
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F(x) 



IX 



tix) 



¥ , (I 12 . ,, 



+K 2 (V8ix) 



it)+ 
it) }> + c.c. 



(A7) 



With the help of this relation, the expression (E9) for the general solution F(x) can now 
be easily deduced. Using the series representation of the modified Bessel functions ||, it 
can be verified F(x) satisfies the perturbative boundary condition F(0) = — 1. 



APPENDIX B 

We discuss here the determination of the constants Ci,C 2 in the infrared limit A — > 0, 
and obtain the asymptotic forms of F(x) and G(x) for large values of x. 

To this end we substitute (p5|) and (|26|) back into the original integral equations (j22"D and 
( p3|) . When the expression fl2~9|) for is used in these equations, we obtain two relations 
which must be satisfied identically for all values of x. The constants C\ and C 2 are fixed by 
the consistency conditions which arise from these relations at the points x = ^P- and x = a, 
respectively. Using the properties of the modified Bessel functions , it is straightforward to 
show in the limit A — ► 0, that these conditions demand C\ to vanish. Then C 2 is determined 
by the system: 

x^d_ ( G{x) \ { 
2 dx \ x 2 I 

\ / x=a 

-2te[-*r) ~ H + 8^ + ^ (B1) 

\ / x=a 

This set of equations fixes the complex constant C 2 uniquely, since the discriminant D of 
the system formed by C 2 and C 2 is non- vanishing. To evaluate it, we use the expression ( |2T)| ) 
which gives the general solution F(x) in terms of the modified Bessel functions, together 
with the fact that C% = 0. Futhermore, we express G(x) [eq. (p5|) 1 in terms of F(x) and 
its derivatives, as well as H [eq. (|24"D l in terms of an integral over F(x). These factors 
make the explicit expression of the discriminant have a rather complicated form involving 
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the modified Bessel functions. Using the relations satisfied by these functions, we find after 
a straightforward calculation that: 



D = -ia [ 2a + -j- ) K 2 {(3)K 2 {(3*) - ia 2 {a - 11) K' 2 {(3)K' 2 {(3*) -\ 



+i 



4 

a /33 3 



( Ilia + ^-a + Aia 2 -11a 2 ) K' 2 {P)K 2 {(3*) +c.c. 



(B2) 



where (3 = y/8ia and K' 2 (z) = -£K 2 (z). 

Since C\ vanishes, we see from equation (^) that the complementary solution decreases 
exponentially for large values of x. Then, the asymptotic form of F(x) is determined by 
F(x), which in the infrared limit A — ^ 0, becomes: 

rx df poo (j£ ~) 

K 2 (V8ix) j —I 2 (V8it) + I 2 (V8ix) J —K 2 (V8it) \ + c.c. (B3) 

In order to evaluate F(x) for large values of x, we use the asymptotic expansions of the 
modified Bessel functions [[|: 



K M = M » i] (B4) 



V2nz 

Then, the leading contributions arising from the integrations in (|B3|) can be calculated 
explicitly. In this way, we find that the asymptotic behavior of F(x), for large values of x, 
is given by: 

F(x) - [x » 1] (B5) 



With the help of this result, we obtain from (|25|) the asymptotic form of G{x) for large 
values of x: 

3Q 

G(x) - -- [x » 1] (B6) 

Since x = these expressions describe equivalently the behavior of our solutions for 
small values of the momenta. 
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APPENDIX C 



Here we present a very simple model which illustrates a possible mechanism for the 
cancellation of the infrared divergences, outside the ladder approximation. The model has 
some similarities with the magnetic sector of the thermal Yang-Mills theory. We consider 
the following expression for the pressure: 

p = p Q + cTY,f ^^nfepo) (ci) 

n J P 2 + Po 

The perturbative series of Il(p — > 0,p = 0) is given by: 

fl(p) = -ag 2 P T + b(g 2 T) 2 + ... (C2) 

where a, b and c are positive constants. We see that a thermal mass of order g 2 T is generated, 
a possibility which might occur also in the case of the magnetic mass. Our model assumes, 
for definiteness, that the contributions to H(p) exponentiate. Then, flU2p can be written as 
follows: 

fi(p) = ^/ (e-? - l) (G3) 

where a = — q 2 . 

Proceeding in parallel with previous cases, we obtain that the expression describing the 
dominant infrared behavior of the pressure has the form: 

P ^2n^T /Ve"^ dp (C4) 
b J\ 

Using this form, we find that the perturbative expansion of the pressure contains power 
infrared divergences given by: 



b 



a 3 , fT\ a 4 T a 5 /T\ 2 
' 3! lU UJ + ¥ A ~ 5!2 IaJ +••• 



(C5) 



Note that this behavior is similar to the one shown by the perturbative thermal Yang- 
Mills theory. 
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However, the full expression ( |U4|) is well behaved in the infrared limit A — > 0. Indeed, 
introducing the variable x = — , we find in this limit that P becomes: 



V 



P 2tt^TV / — e~ x (C6) 



6 ./a x 4 

This can be evaluated in terms of the exponential-integral function 



P ^ !L^r 4 

"3 4 



a 3 Ei(-a) + (2 - a + a 2 )e~ a (C7) 



Therefore, in our model the complete expression giving the pressure in the infrared 
domain in a well behaved function of a. 
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FIGURES 

FIG. 1. The Schwinger-Dyson equation for the scalar self-energy function in the ladder approx- 
imation. 

FIG. 2. (a) A general contribution of the thermal scalar field to the pressure in the ladder 
approximation and (b) its graphical representation as given by equation 

FIG. 3. Graphical representation of equation (|l(^). Wavy lines denote bare gluons and dashed 
lines represent bare ghost particles. Combinatorial factors are shown explicitly in the diagrams. 

FIG. 4. (a) The Schwinger-Dyson equation for the gluon polarization tensor in the ladder 
approximation and (b) the corresponding equation for the ghost two-point function. 



19 



